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ABSTRACT 

We present a spectrum of the symbiotic star V1016 Cyg observed with the 3.6 m Canada-France- 
Hawaii Telescope, in order to iUustrate a method to measure the covering factor of the neutral scattering 
region around the giant component with respect to the hot emission region around the white dwarf 
component. In the spectrum, we find broad wings around Ha and a broad emission feature around 
6545 A that is blended with the [N II] A 6548 line. These two features are proposed to be formed by 
Raman scattering by atomic hydrogen, where the incident radiation is proposed to be UV continuum 
radiation around Ly/3 in the former case and He II A 1025 emission line arising from n = 6 — ^ n = 2 
transitions for the latter feature. We remove the Ha wings by a template Raman scattering wing profile 
and subtract the [N II] A 6548 line using the 3 times stronger [N II] A 6583 feature in order to isolate 
the He II Raman scattered 6545 A line. We obtain the flux ratio -^6545 /-Feseo = 0.24 of the He II A 6560 
emission line and the 6545 A feature for V1016 Cyg. Under the assumption that the He II emission 
from this object is isotropic, this ratio is converted to the ratio <I'6545/'I'i025 — 0.17 of the number of 
the incident photons and that of the scattered photons. This implies that the scattering region with H I 
column density Nhi > 10^" cm~^ covers 17 per cent of the emission region. By combining the presumed 
binary period ^ 100 yrs of this system we infer that a significant fraction of the slow stellar wind from 
the Mira component is ionized and that the scattering region around the Mira extends a few tens of 
AU, which is closely associated with the mass loss process of the Mira component. It is argued that the 
Raman scattered He II 6545 line is an important and useful tool to investigate the mass loss process 
occurring in the late stage of stellar evolution. 

variable star (Mira) — binaries : symbiotic — stars : individual (V1016 



Subject headings: scattering 
Cyg) 

1. INTRODUCTION 

V1016 Cygni is known to be a symbiotic star consist- 
ing of a Mira variable and a white dwarf. Mira variables 
usually exhibit a heavy mass loss with a typical value 
M = 10^^ M0 yr~^ in the form of a slow stellar wind 
with a wind speed ~ 10 km s~^. V1016 Cyg showed an 
outburst in 1965, which appears to be a slow nova-like 
eruption, indicative of a mass transfer process between the 
Mira and the white dwarf (Fitzgerald et al. 1966). The 
outburst is proposed to be thermonuclear runaway on the 
surface of the white dwarf component. 

In symbiotic stars, strong UV radiation, which is pro- 
duced around the mass accreting white dwarf, ionizes the 
surrounding nebula giving rise to prominent emission lines 
(e.g. Kenyon 1986). Furthermore, a significant part of a 
slow stellar wind from the Mira component will be ion- 
ized, whereas a predominantly neutral region may reside 
in the vicinity of the giant star. Therefore, useful infor- 
mation about binary interaction may be obtained from an 
estimate of the extent of the neutral region around the 
giant. 

About a half of the symbiotic stars including V1016 Cyg 



are known to exhibit broad emission features around 6825 
A and 7088 A, which Schmid (1989) identified as the Ra- 
man scattered O VI AA 1032, 1038 by atomic hydrogen. 
The Raman scattering process starts with a scattering hy- 
drogen atom initially in the Is state and an incident UV 
photon redward of Ly/3 and ends up with the hydrogen 
atom in the 2s state and an optical photon redward of Ha. 
The operation of this Raman scattering process requires 
the co-existence of a hot ionized O VI emission nebula and 
a cold neutral scattering region, which is plausibly satis- 
fied in the case of symbiotic stars. Simultaneous detections 
of the features in the far UV and optical regions support 
strongly the Raman scattering nature of these emission 
features (Birriel, Espey & Schulte-Ladbeck 1998, 2000). 

High resolution spectroscopy and spectropolarimetry 
show that these Raman scattered O VI features exhibit 
multiple peak structures and strong polarization often ac- 
companied by polarization flip in the red part (Schmid & 
Schild 1994, Harries & Howarth 1996, 2000). Whereas Ra- 
man scattering of O VI 1032, 1038 operates in the presence 
of neutral H I scattering regions with H I column density 
Nhi ~ 10^^ cm~^, UV incident radiation much closer to 
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Lyman series is scattered in a region with much less Nhi- 
Examples of these processes include the formation of broad 
Hq wings and He II Raman scattered lines. 

Nussbaumer, Schmid & Vogel (1989) advanced an idea 
that Ha wings can be made through Raman scattering 
of Ly/3, which was applied in a quantitative way by Lee & 
Hyung(2000) in their analysis of the Ha wings of the plan- 
etary nebula IC 4997. Many symbiotic stars also exhibit 
similar broad Ha wings that are fitted very well with the 
template wing profiles obtained from Raman scattering of 
Ly/3 (Lee 2000, see also Ivison, Bode, & Meaburn 1994, 
van Winckel, Duerbeck, & Schwarz 1993). 

With their recent extensive spectroscopy of young plan- 
etary nebulae and post AGB stars, Arrieta & Torres- 
Peimbert (2003) found that Ha lines exhibit broad wings 
in these objects. They showed that the wings are well fit- 
ted by a simple formula proportional to AA~^, giving sup- 
port for the Raman scattering origin of Ha wings. Broad 
wings around Ha are also reported in many post AGB ob- 
jects (e.g. Balick 1989, Sclvclh & Bonafacio 2000, van de 
Steene, Wood & van Hoof 2000, Arrieta & Torres- Peimbert 
2000). 

A He II ion being a single electron atom. He II emis- 
sion lines arising from the transitions between energy lev- 
els with even principal quantum numbers have very close 
wavelengths to those for the H I resonance transitions. 
Van Groningen (1993) found a broad feature around 4851 
A in the spectrum of the symbiotic nova RR Telescopii, 
which he proposed to be formed through Raman scatter- 
ing of He II emission lines by atomic hydrogen. Raman 
scattered He II may be observed in broader classes of ob- 
jects than symbiotic stars. In particular, in the planetary 
nebula NGC 7027, Pequignot et al. (1997) reported the 
existence of He II Raman-scattered feature around 4851 
A blueward of H/3. He II Raman-scattered 6545 line was 
also found in the symbiotic stars RR Tel and He2-106 (Lee, 
Kang & Byun 2001). 

In this paper, we present a spectrum of the symbiotic 
star V1016 Cyg observed with the 3.6 m Canada-France- 
Hawaii Telescope, from which we isolate the He II Ra- 
man scattered 6545 A feature. We calculate the ratio of 
the number fluxes of incident and scattered radiation in a 
straightforward way, which will provide a direct measure 
of the covering factor of the neutral scattering region. 

2. DATA AND ANALYSIS 

We observed the symbiotic star V1016 Cyg on the night 
May 20, 2002 with the coude spectrograph 'Gecko' at the 
3.6 m Canada-Prance-Hawaii Telescope at Mauna Kea. A 
2kx 4.5k EEVl CCD chip was used for the detector, and 
the order sorting was made with an interference filter. The 
spectral resolution was R — 120, 000. Using IRAF pack- 
ages, standard procedures have been followed in order to 
reduce the data. 

In Fig. 1, we show the spectrum of V1016 Cyg. We show 
fluxes in the vertical axis with arbitrary unit, because in 
this work the only necessary quantities are relative ratios. 
The vertical arrows in the figure mark He II A 6527, the 
Raman scattered He II 6545 A, and He II A 6560. We blow 
up the data by a factor of 50 in order to see clearly the 
weak He II A 6527 and the 6545 feature represented by the 
solid thick line. 



The emission line He II A 6527 arises from a transition 
between n = 14 and n = 5 states of He II. Raman scat- 
tered He II 6545 A feature is formed when a He II A 1025 
line photon (n = 6 ^ n = 2) is incident upon a hydrogen 
atom in the ground state followed by a de-excitation of 
the scattering hydrogen atom into 2s state re-emitting an 
optical photon with A = 6545 A. Because of the energy 
conservation, the wavelength Ao of the Raman scattered 
radiation is associated with that of the incident UV radi- 
ation by 

(1) 



Ao — (Aj ^ - Xi^a) 



where XLya is the wavelength of hydrogen Lya. The varia- 
tion AAi in the incident radiation results in corresponding 
variation of AAo given by 

AA„ A„ AA, 



A,; A,; 



which leads to the broadening of the Raman scattered ra- 
diation by a factor of 6.4. 

He II emission lines from 2n — > 2 transitions have wave- 
lengths that differ from those of H I Lyman series lines 
by 1/2000 due to the small difference in reduced masses. 
This results in Raman-scattered radiation with wavelength 
differences from Balmer lines 



AAo ^ 5.9 



n2(n2 - 1) 



(n2 - 4)2 



A, 



(3) 



(Lee et al. 2001). 

The scattering cross section for this process can be com- 
puted using the Kramers-Heisenberg formula, which yields 
fj ~ 10^20 cm-2 g a,l. 2001, Nussbaumer et al. 

1989). This cross section is two orders of magnitude larger 
than the corresponding values for the Raman scattering of 
O VI 1032, 1038. Fig. 2 shows the energy level diagrams 
of hydrogen and He II, where emission lines of He II and 
the Raman scattering processes are illustrated. 

In this spectrum, there appear very broad wings around 
Ha, which are excellently fitted in the far wing parts by 
AA^2 = (^ _ A_f/Q)~2 profile represented by a dotted line 
in Fig. 1, Aija being the difference of the wavelength from 
that of Ha. This profile is proposed to be originated from 
Raman scattering of fiat incident radiation around Ly/3 
by Lee (2000) and Lee & Hyung (2000). We prepare the 
template profiles of Raman scattering wings that are ex- 
pected to be formed in neutral scattering regions with H I 



column densities Nhi 



10 



19 



10 cm irradiated by 



a flat continuum radiation around Ly/3, as was done by 
Lee & Hyung (2000). These template wings are shown in 
Fig. 3. The exact choice of the H I column density does 
not alter the main result of this work, because it affects 
the shape near the Ha center and we arc only interested in 
the relative strength ratio of He II emission and Raman- 
scattered lines that lie in the far blue wing part of Ha. 
We select the template wing profile for Nhi = 10^^ cm"^^ 
which is also capable of producing a Raman scattered He II 
6545 feature. 

After subtraction of the Ha wings, we perform Gaus- 
sian fittings with a functional form /(A) = /oexp[— (A — 
Ao)V^A2] to He II emission lines at 6527 A and 6560 A, 
and the Ha, using a least square method. Satisfactory 
fits are obtained by a single Gaussian for He II emission 
lines at 6527 A and 6560 A and Ha. The line widths 
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are AA = 0.96 A, 0.98 A for He II 6527 and 6560 respec- 
tively, and for Ha we obtain AA = l.oA. The full widths 
at half maximum of these emission lines are 73 km s~^, 
74 km s~"^ and 78 km s~^, respectively. The peak values 
arc /o = 5.2 x 10^ for He II A 6560 and /o = 1.8 x lO'"^ for 
He II A 6527. Fig. 4 shows our results of Gaussian fittings 
for H I and He II emission lines. 

In order to separate the Raman scattered He II 6545 
feature from the spectrum, we make use of the fact that 
the [N II] A 6548 profile slioTilcl be identical with that of 
the [N II] A 6584 feature, which is 3 times stronger due 
to the atomic structure (e.g. Osterbrock 1989, Storey & 
Zeippen 2000). We subtracted the flux of the [N II] A 6584 
multiplied by 1 /3 (represented by red lines) from the flux 
around the 6545 and [N II] A 6548 features (represented 
by the black lines) to isolate the He II Raman-scattered 
feature. 

We show the result in Fig. 5, where the Raman scat- 
tered He II 6545 feature is seen clearly in solid thick line. 
The dotted line is a Gaussian fit to the feature, where the 
width is AA = 6.2 A and the peak value is /o = 1.0 x 10^. 
This width is simply chosen from the mean value of the 
widths of the He II emission lines multiplied by a factor 
6.4, which is the enhancement factor of the scattered line 
width attributed to the scattering incoherency. 

3. CALCULATION 

3.1. Strength of Raman Scattered He II 654-5 Line 

The He II A 1025 flux of V1016 Cyg can only be in- 
directly inferred from the fluxes of He II emission lines 
at 6527 A and 6560 A. From the recombination theory, 
relative strengths of He II emission lines can be calcu- 
lated as a function of temperatures and electron densi- 
ties (e.g. Storey & Hummer 1995, Brocklehurst 1971). 
According to Gurzadyan (1997) the calculated fluxes of 
He II A 6560 and He II A 1025 lines are given by i^eseo = 
0.135F4686, FiQ2b = 0.618i^4686, where F4686 is the flux 
of He II 4686. In this work, it is assumed that the 
temperature is T = 20000 K and the electron density 
rif. = 10* cm^"' for illustrative purpose. If we consider 
the case B recombination, the ratio of the number fluxes 
of the He II A 6560 and He II A 1025 line photons generated 
from recombination is 

*io25/$6560 = (0.6 1 8 /0.135)(1025/ 65 60) = 0.7 1 5. (4) 

This implies that the number of He II A 1025 line photons 
is almost 70 per cent that of He II A 6560 line photons, 
which does not vary sensitively with the physical condi- 
tions assumed above. Because the Raman He II 6545 is 
located close to He II A 6560, their flux ratio directly in- 
dicates the scattering efBciency. 

Under the assumption that the He II emission from this 
object is isotropic, this ratio is converted to the ratio of the 
number of the incident photons and that of the scattered 
photons 



^6545/^1025 



1.1 X 10^ 6545 $1025 

5.2 X 104 1025 $ 



6560 



0.17. (5) 



Birriel et al. (2000) used the Hopkins Ultraviolet Telescope 
(HUT) and various ground based telescopes to measure the 
strengths of the resonance doublet O VI 1032, 1038 and 
their Raman scattered features at 6825 Aand 7088 A, in 



order to obtain the Raman scattering efflcicncies in a num- 
ber of symbiotic stars. They conclude that the scattering 
efficiency for V1016 Cyg is about 50 %, which is larger 
than our own estimate. 



3.2. Neutral Scattering Region and Mass Loss Process of 

the Mira 

3.2.1. Covering Factor of the Neutral Scattering Region 

In this subsection we consider the neutral scattering re- 
gion in V1016 Cyg which is responsible for the formation 
of the 6545 feature. The giant component of V1016 Cyg is 
believed to be a Mira variable, and the white dwarf com- 
ponent showed outburst in 1965. A Mira variable usually 
exhibits a mass loss at a rate M ^ 10^^ _ 10^^ yr^-"^ 
in the form of a slow stellar wind with a typical wind veloc- 
ity t; ~ 10 km s~^ (e.g. Vassiliadis & Wood 1993, Schild 
1989). The neutral scattering region is plausibly formed 
around the giant component and closely related with the 
mass loss process of the giant. 

The mass loss rate of a Mira variable may not be a con- 
tinuoTis process and can be episodic as it evolves in the 
asymptotic giant branch, which was proposed by Zijlstra 
et al. (1992). However, for simplicity, we assume that 
the mass loss rate is constant and the mass loss process is 
spherically symmetric. In a binary system like V1016 Cyg, 
this assumption may be quite unrealistic because the ma- 
terial distribution will be severely affected by the presence 
of the white dwarf companion and the orbital motion. As- 
suming a wind consisting of only hydrogen we may write 
the mass loss rate 



M = 47rr^mpn(r)w(r). 



(6) 



The slow stellar wind velocity is almost constant in a 

quite extended region, and therefore the density is ex- 
pected to decrease as r^^. If we assume that this spher- 
ically symmetric region consists of complete neutral hy- 
drogen, any sight line originating from the emission re- 
gion around the white dwarf will have the Raman scatter- 
ing optical depth exceeding unity for a mass loss rate of 
M = 10~^ Mq yr~^. Therefore, we may expect that the 
wind material should be ionized significantly. The ioniza- 
tion front is generally sharply defined on the order of the 
mean free path of Lyman limit photons, which is quite 
small compared with the characteristic scale of the scat- 
tering region. 

Schmid (1995) considered the distribution of neutral hy- 
drogen around the giant component when he investigate 
Rayleigh and Raman scattering processes. In his com- 
putation, he introduced a dimensionless parameter Xh+ 
defined by 



- 



aBiH.T„)aj, 



-DL 



H 



M 



(7) 



Here, aB{H,Te) is the recombination coeflicient of hydro- 
gen in a region with electron temperature T^, D is the 
binary separation, Lh is the number of ionizing photons 
from the hot white dwarf, jj, is the mean molecular weight 
and ttH is the rclatvie hydrogen particle density. 

For given 9 the distance s = uL to the ionization front 
from the He II emission region is computed by considering 
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the integral given by 
Jo 



X dx 



(a;2 -2a; cos (9 + 1)2 

1 \1 { T U — COS 6 „ IT 

- tan"^ — e+- 

2 V sin6' 2 



sin'^ 



1 . „ / -1 w — cos( 

+ - sin 2 tan : — — 

4 V sin ^ 



6* 



(8) 



The ionization boundary can be located when Xfj+ ex- 
ceeds f{9) defined by 

f{0)= lim f{9,s) 



S^OO 

1 

sin^ 



ir-e sm20 
+ —r- 



(9) 



Fig. 6 shows the graph of f{9)- As is pointed by Schmid 

(1995), when Xh+ > 7r/4 = f(0 = 7r/2), most of the giant 
wind is ionized leaving a conc-likc neutral region around 
the giant. 

It is expected that the sight line from the He II emission 
region will be characterized by a large H I column density 
Nhi > 10^^ cm~'^, at which the Raman scattering opti- 
cal depth becomes unity. Therefore, wc may safely regard 
that He II A 1025 line photons will be Raman-scattered 
with almost same efficiency. Lee & Lee (1997) performed 
detailed Monte Carlo calculations to show that the effi- 
ciency is about 0.6, which means that about 60 per cent of 
He II A 1025 photons that enter a neutral region with the 
Raman scattering optical depth exceeding unity will come 
out as Raman-scattered optical photons. The remaining 
40 per cent of UV photons are reflected near the surface 
by Raylcigh scattering. Since 17 per cent of He II A 1025 
photons are converted into Raman scattered He II 6545 
photons, it follows that the covering factor Cg of the scat- 
tering region is Cs ~ 0.17/0.6 = 0.28. Noting that 

1^ (10) 



Afl 

477 



we obtain 6s = 64°, where Ail is the solid angle of the 
scattering region subtended with respect to the He II emis- 
sion region. Substituting 9s = 64° into Eq. (9), we obtain 



Xh+ = f{0s) ^ 1.9. 



(11) 



This estimate is similar to the value Xh+ 2.3 — 3.7 pro- 
posed by Miirset et al. (1991), but quite large compared 
with the value Xh+ ~ 0.4 obtained by Birriel et al. (2000). 

3.2.2. Mass Loss Rate of the Mira 

Despite intensive research over several decades, the or- 
bital period of V1016 Cyg is only poorly known. The 
distance to V1016 Cyg is also highly uncertain with esti- 
mates ranging from 1 kpc to 10 kpc (Watson et al. 2000 
and references therein). From a recent analysis using opti- 
cal images obtained with the Hubble Space Telescope and 
radio images from the Vary Large Array (VLA), Brock- 
sopp et al. (2002) suggested a projected binary separation 
of 84 ± 2 AU with the assumption of a distance 2 kpc. Us- 
ing this separation, they proposed that the binary period 
will be longer than 500 years. Previous works by other 
researchers provide smaller binary separations and shorter 
binary periods. For example, Schmid & Schild (2002) per- 
formed spectropolarimetry to investigate the variation of 



the polarization direction which is expected to vary due to 
the orbital motion. Based on their polarimetric method, 
they proposed that the binary period may exceed 100 
years. 

Here, wc assume that the distance to V1016 Cyg is 2 kpc 
as Brocksopp et al. (2002) did. With the H I distribution 
given in Eq. (6), the H I column density for a sight line 
into the neutral region is nearly inversely proportional to 
the impact parameter p, 

- (I) (t) 

where i?* , n, is the radius of the Mira and the hydrogen 
density near the surface of the Mira. Assuming the mass 
loss given in Eq. (6) we have 



n*ii* ( — 1 = 2.5 X 10 
P 



21 —2 

cm 



Mr. 



1 AU 
P 



'10 km s-i 



, (13) 

IO-^Mq yr-1 ' " " ' 

The impact parameter p for the sight line that grazes the 
ionization boundary with 9 = 6s, we have p = Dsm9s = 
0.9D, where D is the binary separation of V1016 Cyg. The 
Raman scattering optical depth for this sight line exceeds 
unity, so that Nhi > 0.9 x 10^^ cm"^. Therefore, we have 

If we adopt the binary separation D = 80 AU obtained 

by Brocksopp et al. (2002), and assume the stellar wind 
velocity about — 10 km s^^ near the Mira surface, we 
obtain M > 3 x 10~^ Mq yr~^. However, this estimate is 
highly dependent on the choice of the distance to V1016 
Cyg, which is currently quite uncertain. 

4. DISCUSSION 

In this paper, we introduced a method to obtain a cov- 
ering factor of the neutral scattering region with respect to 
the emission region formed around the white dwarf com- 
ponent in a symbiotic star. Since the spacing between the 
Raman scattered 6545 feature and He II A 6560 emission 
line is quite small, a single exposure is needed to obtain 
both features and the effects including interstellar extinc- 
tion can be avoided as far as we are concerned with the 
flux ratio of both features. Furthermore, the Raman scat- 
tered He II 6545 feature can be isolated with relative ease 
from [N II] A 6548 feature using always 3 times stronger 
[N II] A 6584 line, which makes the current method quite 
robust. With this ratio, we may obtain important infor- 
mation regarding the ionization radiation and the mass 
loss process in symbiotic stars. 

Both Ha broad wings and the 6545 feature have Raman 
scattering origin, where continuum around Ly/3 is respon- 
sible for the former feature and He II A 1025 emission line 
is for the latter. Therefore, the strength ratio of the 6545 
feature relative to that of Ha wings is equal to the flux 
ratio of the He II A 1025 line to that of the continuum 
around 1025 A. Since the Ha wings are approximated by 
1.4 X 10'(A — XHa)^^, and the 6545 feature is fitted with 

the Gaussian 1.12 x iQ^e-^'-^-^^'^^'^y'-'^^^^^ the equivalent 
width of the 6545 feature 

We545 = ( / e-(^/^ ^^'dx) A = 21 A (15) 
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Noting that the wavelength scale is enlarged by a factor of 
6.4 in the optical region due to the incoherence of Raman 
scattering, we obtain the equivalent width of He II A 1025 
is EWHeiiw25 = 3.3 A. A direct measurement of this 
quantity is severely hindered by the interstellar medium. 
However, we may infer it by measuring the equivalent 
width of nearby He II 1080 Une using FUSE. 

Espey et al. (1995) used the data obtained with the 
HUT to estimate that 25-50 per cent of O VI 1032, 1038 
doublet photons interact in the neutral scattering region 
in the symbiotic star RR Tel. Because the distance and 
the binary orbital period of RR Tel are also very uncer- 
tain, a direct comparison with the results for V1016 Cyg is 
not possible. However, it appears that a very similar scat- 
tering geometry is shared by these two symbiotic novae. 
Even though He II A 1025 has two orders of magnitude 
larger scattering cross section than O VI 1032, 1038 lines, 
it may be that they share almost same neutral scatter- 
ing region due to strong H-ionizing UV radiation from the 
white dwarf component. 

Schmid et al. (1999) proposed that the strength of 
Raman scattered O VI 6827, 7088 may change with the 
binary orbital phase because the scattering process is a 
dipole process where forward and backward scattering is 
more favored than orthogonal scattering. The same argu- 
ment also applies to the Raman scattered He II 6545 lines, 
where the flux of He II 6545 feature is expected to be 
stronger in the conjunction phases than in the quadrature 
phases. In order to test this expectation, it requires a very 
long term monitoring. More information may be obtained 
from spectropolarimetry as has been done for Raman scat- 
tered O VI lines (e.g. Harries & Howarth 1996, Schmid & 
Schild 1996). Binary orbit parameters have been extracted 
in some symbiotic stars by noting the variation of the po- 
sition angle. 

There may be additional scattering regions that are not 
considered in this work. For example, spectropolarimetry 
around Raman scattered O VI 6827 in many symbiotic 
stars shows that the polarization direction in the reddest 



part with v > +80 km is flipped by an amount of al- 
most 90° with respect to the strongly polarized blue and 
main part (e.g. Harries & Howarth 1996). Schmid (1996) 
proposed that this flip is attributed to the receding part of 
the slow stellar wind around the giant component. How- 
ever, the typical wind speed is ~ 10 km s~^, which is much 
smaller than the observed values. 

Lee & Park (1998) proposed that the profile of O VI 
6827 is mainly attributed to the kinematics of the emis- 
sion region formed around an accretion disk. The emission 
region may extend in sub-AU scale with the Keplerian 
speed ~ 100 km s~^ in the orbital plane. In this picture, 
the flipped red part can be explained by proposing the 
existence of additional neutral or clumpy regions that are 
moving away in the direction normal to the orbital plane. 
This raises an interesting possibility that this additional 
component is related with the bipolar morphologies exhib- 
ited in many symbiotic stars. However, this additionally 
scattered component appears weak in O VI 6827 features, 
and maybe this does not contribute much to the He II 
6545 feature. Future spectropolarimetry using big tele- 
scopes may shed more light about this. 

We conclude that the Raman scattered He II 6545 fea- 
ture is extremely useful to investigate the mass loss process 
of a giant star and the size of the neutral scattering region, 
which will put very tight constraints on the distance and 
binary orbital period of symbiotic stars. 
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Raman Scattered He II 6545 in V1016 Cyg 
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Fig. 1. The spectrum around Ha of the symbiotic Mira V1016 Cyg obtained with the 3.6 m Canada- France-Hawaii Telescope. The 
vertical arrows mark He II A 6527, the Raman scattered He II 6545 A, and He II A 6560. The solid thick line is a blow-up of the solid thin 
line by a factor of 50, where the weak He II A 6527 and the 6545 feature are clearly illustrated. The flux is shown in arbitrary units because 
physically meaningful quantities are relative fluxes in the current work. Note that there appear very broad wings around Ha, which are 
excellently fitted by AA~^ = (A — XHct)~'^ profile represented by the dotted line, Xhoi being the difference of the wavelength from that of Ha. 
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Fig. 2. — Atomic level structures of hydrogen and He II. He II emission lines at 6527 A and 6560 A arise from n = 14 — ► n = 5 transitions 
and n = 6 — > n = 4 transitions, respectively. He II A 1025 (n = 6 — > n = 2) emission line has a slightly higher energy than Ly/3 1025. When a 
He II A 1025 photon is incident upon a hydrogen atom, it may de-excite to 2s state by re-emitting an optical photon blueward of Ha which 
appears at 6545 A. 
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Fig. 3. - Template wing profiles around Ha that are formed from Raman scattering of a flat continuum around Ly/3 in a neutral scattering 
region with H I column density Nhi- The wing profiles were prepared using a Monte Carlo technique, of which the detailed procedure is 
explained in Lee & Hyung (2000). 
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Fig. 4. — Gaussian fits using a functional form /(A) = /o cxp [-(A - A())2/AA2] to He II A 6527, He II A 6560 and Ha after subtracting 
tfie wide Hq wings excellently approximated by AA~^. Because Hq is too strong and He II A 6527 is too weak, we present the results in two 
panels. In the upper panel, we show the Gaussian fits to He II A 6527 and Ha by dotted lines, where all the emission lines are nicely fitted 
with a single Gaussian. The peak values are fo = 1.7 X 10^ for He II A 6527 and fo = 4.8 x 10* for He II A 6560. The lower panel shows the 
detailed result around He II A 6560. 
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Fig. 5. — The Raman scattered He II 6545 isolated by deblending from the [N II] A 6548 emission line. Because forbidden lines arc optically 
thin and [N II] A 6548 and [N II] A 6584 originate from the same excited level, the profiles arc exactly same and differs only in strength by 
a factor of 3. The dashed line shows the [N II] A 6584 emission line divided by 3 and translated to the location of [N II] A 6548. The dotted 
line is a single Gaussian fit to the Raman scattered He II 6545 feature. The peak value /o = 1.0 x 10^ and the width is 6.4 times that of the 
He II A 6527 and He II A 6560, which is exactly expected from the incoherency of Raman scattering. 
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Fig. 6. — The scattering geometry of V1016 Cyg and the ionization of the stellar wind from the Mira. Due to the strong UV radiation 
coming from the mass accreting white dwarf, a significant part of the slow stellar wind will be ionized. The thick curve shown in the upper 
panel represents the ioiiizatiou front dividing the neutral region around the giant and ionized region formed in the part where the white dwarf 
resides. The low panel shows the behavior of the function f(9) defined in the text. 



